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Abstract 

The core-core structure factor of dense star polymer solutions in a good sol- 
vent is shown theoretically to exhibit an unusual behaviour above the overlap 
concentration. Unlike usual liquids, these solutions display a structure fac- 
tor whose first peak decreases by increasing density while the second peak 
grows. The scenario repeats itself with the subsequent peaks as the density 
is further enhanced. For low enough arm numbers / (/ < 32), various differ- 
ent considerations lead to the conclusion that the system remains fluid at all 
concentrations. 

PACS: 82.70Dd, 61.20Ja, 61.20.Gy, 61.25Hq 
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I. INTRODUCTION 



Star polymers are macromolecular entities consisting of / polymeric arms chemically 
attached to a common centre. At the limit where the number of monomers per arm (degree 
of polymerisation) is large, the size of the central core of the star is negligible compared 
with the overall radius of the macromolecular aggregate and can be ignored at a first approx- 
imation. Scaling theory has provided insight into the structure and conformation of a single 
star and the way that these properties are influenced by the arm number (functionality) /, 
N and solvent quality From both theoretical and experimental points of view, it is 

even more interesting to describe the properties of dense star polymer solutions by means 
of an effective pair potential between star centres. Experimentally, this is important as 
interstar correlations are ubiquitous in understanding neutron scattering spectra of dense 
star polymer solutions. Moreover, such a description of stars establishes a bridge between 
polymer and colloidal science: indeed, due to their peculiar construction, star polymers can 
be viewed as hybrids between polymers and colloidal particles. Once such a description 
has been established, it is then possible to look at a star polymer solution as an effective 
one-component system of point particles (the star centres) and employ the known machinery 
from liquid state theory and/or computer simulations to study its properties. 

By way of direct comparison with experimental data, it was found recently that star 
polymers in a good solvent can be described by an effective pair potential which is logarith- 
mic for short distances and crosses over to a Yukawa form for larger interstar separations 
0]. This is a new type of interaction in the sense that (i) it has an ultra-soft logarithmic 
core whose hardness depends in fact on the functionality / in a way which will be explained 
shortly and (ii) it features a crossover from one functional form to another at some length 
scale a which is of the order of the corona diameter of the star. Systems which are de- 
scribed by simple, spherically symmetric interactions have been studied extensively in the 
last thirty years, since the development of powerful computers made integral equation theo- 
ries and simulations computationally tractable [^^. For most of the commonly considered 
interactions (power laws, Yukawa, square- well or square-shoulder potentials etc.) the liquid 
structure factor S{k) features a single characteristic length which is basically set by the 
density. Moreover, the correlations in structure grow with increasing density p until, at 
some temperature-dependent value, the system crystallises. We call the liquids for which 
such a scenario materialises usual or normal. The purpose of this paper is to show that star 
polymer solutions are unusual, in the sense that S{k) or its real-space counterpart, the radial 
distribution function g{r) have quite unexpected behaviour above the overlap concentration 
of the solution. The latter is defined as the polymer concentration at which the stars start 
overlapping within their corona. There exist two competing length scales which manifest 
themselves in the form of a structure factor whose peaks do not grow simultaneously in height 
as p increases, but rather lower order peaks grow higher until they completely dominate and 
the original main peak disappears. These two lengths are the average interparticle distance 
a oc p~^/^ and the corona diameter a of the stars. Moreover, for low enough functionality, 
freezing does not take place, i.e., the system remains fluid at all densities. 

The rest of the paper is organised as follows: in section 2 we present the model pair 
potential and briefly discuss its properties. In section 3 we present results from computer 
simulations and integral equation theories regarding the structure of the system for a very 
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wide range of densities demonstrating the unusual features and we discuss their origin. In 
section 4 we present a simple 'Hard-Sphere mapping' argument to establish the limit of / 
beyond which the system does not crystallise. In section 5 we discuss the connection with 
experiments. Finally, in section 6 we summarise and conclude. 



II. PAIR POTENTIAL BETWEEN STAR POLYMERS 

Our starting point is the pair potential between two star polymers in a good solvent 
separated by a centre-to-centre distance r, which reads as: 



V{r) ^ { (5/18)/3/2 [- ln(r/a) + (1 + ^ /2)-\ 
ksT \ (5/18)/3/2(l + ^/2)-\a/r) exp[-v^(r - cr)/2a 



r > a]. 



where fcs is Boltzmann's constant and T is the absolute temperature. This is an entropic 
interaction stemming from excluded volume effects between monomers belonging to two 
different stars in a good solvent. It was shown recently that the length scale a has to 
be chosen as twice the distance between the centre of the star and the centre of the largest 
(outermost) blob 0]. The logarithmic form of the interaction for interparticle distances r < a 
results from the arguments of Witten and Pincus p|. For r > a we take an exponential form 
for the interaction with a decay length equal to the largest blob diameter. For details on the 
determination of the overall numerical prefactor of the logarithmic term, we refer the reader 
to Ref. 1^. The amplitude of the Yukawa tail is finally determined by the requirement of 
smoothness of the interaction at r = cr. In what follows we use a as the unit of length and 
introduce a dimensionless 'packing fraction' rj defined as: 

V = Ip^'- (2) 

It has been shown in Ref. that this choice of a pair potential ^ brings about good 
agreement between theory and experiment regarding small-angle neutron scattering (SANS) 
data for a wide range of densities (or polymer volume fraction). Moreover, a similar expo- 
nential 'cutoff' of the logarithmic part of the interaction has also been employed recently 
in a study of the elastic moduli of block copolymer micelles [0. The entropic nature of 



the interaction renders the temperature T irrelevant, as the Boltzmann factor e"^^*^*") is 
temperature-independent. Instead, the arm number / plays now the role of an effective 
inverse temperature. Indeed, in the limit / — >■ the interaction vanishes, whereas at the 
"colloidal limit" / — > cxo we recover the well-known hard sphere (HS) interaction. In fig. 
m we show the pair potential for different values of /. Whereas for large / = 128 and 256 
there is a spectacular change in its behaviour as r crosses a, for low /-values, the potential 
is ultra-soft and its range becomes longer. 

As will be shown in section 4, the ultra-soft character of the interaction at hand, has the 
consequence that, for low enough values of the functionality /, the system never crystallises, 
no matter how large the density or external pressure are. Thus, this system offers us a 
unique possibility to examine the behaviour of the structural functions of the liquid (the 
radial distribution function and/or the structure factor) over a practically unlimited range 
of densities. This is the subject of the following section. 
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III. ANOMALOUS STRUCTURE FACTOR 



First, we summarise some basic notions from the theory of classical liquids and refer the 
reader to Ref. for details. A quantity of central interest for a classical fluid in equilibrium 
is the so-called radial distribution function g{r) and the closely associated pair correlation 
function h{r) = g{r) — 1. If we call p the number density of the liquid, then the quantity 
pg{r) is nothing else but the density profile that develops if we keep a particle fixed at the 
origin. In other words, the quantity g{r) expresses the ordering of the rest of the system 
around a given particle of the liquid. Equivalently, one can look at the structure factor S{k) 
which is simply p times the Fourier transform of the pair correlation function. The peaks of 
S{k) reveal information about the characteristic length scales in the system. Alternatively, 
g{r) is more appropriate if one wishes to determine, e.g., the average coordination number 
of the liquid. 

Both g{r) and S{k) can be measured in a standard Monte Carlo simulation |^,^. We have 
thus performed such simulations for a wide range of densities and, in addition, we have solved 
the Rogers- Young (RY) closure [ill] order to obtain a comparison. After ascertaining that 
the RY gives quite reliable results for all the densities for which simulations were carried 
out, we relied on this closure to calculate the structure of the fluid for very high packings 
(60.00 ^r;^ 10.00), where a simulation would be very expensive since a very large number of 
particles in the simulation box would be necessary in order to obtain reliable results. 

In fig. 1^ we show representative results in order to provide a comparison of the radial 
distribution function g{r) obtained from simulations with the RY result for various different 
values of rj. Apart from small discrepancies for the intermediate value r] = 0.60, it can be 
seen that the agreement is quite good. Hence, the RY closure is a reliable theoretical tool 
for the calculation of the pair structure of the liquid. It should also be noted that for high 
values of r], rj^S.O, the RY closure reduces practically to the Hypernetted Chain (HNC) 
for our system. This is expected since, for high densities, we are dealing with a long-range 
interaction when distances are measures in units of average interparticle distance and the 
HNC is known to be most accurate precisely for long-range potentials. 

The radial distribution function g{r) of the system at hand, shows the typical evolution 
of a normal liquid as the density is increased until we reach the overlap density p*. This 
is the density at which the average interparticle distance a = p^^l"^ becomes equal to the 
length scale a. It corresponds roughly to an overlap packing faction r/* = 0.50. For r/ < rf ^ 
the function gir) displays oscillations with a characteristic length scale a. The heights of 
the peaks of gir^ grow with density; the same behaviour is also observed for the structure 
factor S{k\ of course. 

The situation changes above the overlap concentration, as can be seen in fig. ^(a). Instead 
of having a main peak of g(r^ at r a, we observe that one peak of gir^ always remains at 
slightly outside the logarithmic part at r^^a regardless of the density. The 'main' peak of g{r) 
grows gradually as the density is increased. Indeed, at ?7 = 0.80 g{r) shows just a shoulder 
preceding the peak at r^cr mentioned above and only at ?7 > 1.00 can one distinguish a 
clear peak at r a. Even then, this peak remains lower than the second until rj = 1.30. 
Thereafter, it overpasses the second peak in height but still, the position of the second peak 
of g{r) is not at twice the distance of the first peak, as is the case for normal liquids, because 
the latter is determined by the density and the former by the length scale a. 
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This anomaly in g{r) is reflected, naturally, in the shape and evolution of the structure 
factor S{k), as can be seen in flg. ^(b). As a flrst remark, we flnd that the height of the 
main peak of S{k), which grows until we reach the overlap packing, becomes lower above 
1]*. More unusual features are observed if one looks at the positions and the competition 
between the flrst two peaks. Though the second peak remains at twice the position of 
the flrst up to ?7 = 0.60, at packing rj = 0.80 the position of the flrst peak is practically 
unchanged with respect to that ai rj = 0.60, but the second peak now moves closer to the 
origin and increases in height. This trend persists as r] grows and already at rj = 1.30 the 
second peak has become higher than the flrst, which shows a clear trend of disappearing at 
higher densities, as it will be conflrmed shortly. 

This anomaly in the structure factor can be traced to the crossover of the interaction 
from a Yukawa to a logarithmic form and the softness of the logarithmic potential. Indeed, 
above the overlap concentration the system always tries to maintain one coordination shell 
outside the logarithmic part, where the interaction is weak. At the same time, the softness 
of the logarithmic core allows two things to happen: on the one hand, qualitatively, the 
system remains always fluid which would not have been the case if the core was harder (e.g. 
a HS core, or even the same, logarithmic form for different /, as will be shown in section 
4.) On the other hand, more quantitatively, the softness of the logarithm allows for a rather 
broad flrst peak of the radial distribution function which can accommodate enough particles 
so as to allow for the second peak to occur immediately outside r = a. 

The unusual shape of the structure factor can be understood by means of the existence 
in g{r) of these two different length scales: one length scale a p~^^^, which is manifestly 
density-dependent and one length scale b'^a which is density-independent. Below r]* only 
the flrst length scale appears but above rj* both are present. Let us call kn the position of the 
n-th peak of the structure factor. The flrst peak of S{k) corresponds, roughly, to the length 
scale b, i.e. ki ^ 2-71 /b and the second peak to the length scale b — a, i.e. k2 ~ 27r/(6 — a). 
Indeed, as can be seen from flg. ^(b), ki is practically constant, whereas /c2 decreases with 
density, a feature that can be attributed to the increase of the quantity 6 — a ~ a — p~^^^ 
with growing density. Moreover, the growth of the second peak can be understood since the 
structure that gives rise to it becomes more pronounced with increasing rj, as the flrst peak 
of g{r) takes shape. 

In these terms, we can now make the hypothesis that the flrst peak of the structure 
factor will disappear altogether when the density is such that b — a = aoib = 2a. In 
other words, the structure factor will have a main peak at a position dictated only by the 
density, when the length b is twice the length a in such a way that g{r) has exactly two 
oscillations of wavelength a between r = and r = b. Since 6 ~ a, it turns out that the 
density must be such that a ^ a/2. As can be seen in flg. ^(a) this occurs at the 'magic' 
value r)2 = 3.40, where the subscript denotes the number of oscillations of g{r) in the interval 
[0,6]. Accordingly, the structure factor at r/ = 772 has a strong flrst peak (which, however, 
has evolved from the second peak at lower densities!) located at a position 27r/a ~ 47r/cr. In 
flg. ^(b) we show also the structure factor at 77 = 2.00 in order to show the disappearence of 
what used to be the flrst peak of S{k) and its replacement by the second one. However, we 
emphasize that the fact that at 77 = 772 we have again a structure factor with its main peak 
located at a position k^ax ~ P^^^, does not mean that we are dealing with a normal liquid. 
Indeed, as can be seen from flg. |(b) the structure factor shows some weak substructure 
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which is not present for usual hquids. 

It is now pertinent to ask the question what happens once the density is further increased. 
Does the height of the new main peak grow until the system solidifies or does the above 
scenario repeat itself? We have solved the RY-closure for packing fractions up to 77 = 60.00 
and we find that, in fact, it the second possibility that materialises. As rj grows, the main 
peak is lowered and the one that used to be the third at low densities grows, see fig. H(a). 
The mechanism that brings about this scenario is no other but the development of more 
and more oscillations inside the logarithmic core in the function g{r), as can be seen in fig. 
^(b). In fact, one can repeat the argument about the 'magic' packing fraction r/2 above for 
an arbitrary number of oscillations as follows: the structure factor S{k) will be dominated 
by a single length scale whenever there is an integer number of oscillations m between r = 
and r = h (J. Given the magic packing fraction 772, it is straightforward to show that the 
general magic packing fraction rj^ will be related to ri2 by: 



Using ri2 = 3.40 we obtain rj^ = 11.475, 774 = 27.20 and 775 = 53.125. The structure factor 
for these magic values is shown in fig. It can be seen clearly that for the 7n-th magic 
value, S{k) has a dominant peak located at the position 27rm/a. As the order 7?t, grows, 
the height of the dominant peak also decreases slightly and some substructure develops in 
S{k). However, the main length scale comes from the m oscillations of the radial distribution 
function g{r) in the logarithmic core. 

The above results can be nicely summarised by making a log-log plot of the positions 
of the first few peaks of the structure factor against the packing fraction 7/, as shown in 
fig. 1^. We emphasise here that when we talk about the n-th peak we actually mean 'the 
peak which is the n-th if we extrapolate at low enough densities so that we are in a regime 
where the liquid is normal.' The reason for this distinction is that as 7/ grows the low-order 
maxima of the structure factor disappear, as explained above, and as a result what used to 
be a second-order maximum becomes now first order etc. This is manifested in fig. ^ by the 
fact that the curves representing the positions of the various maxima fc„ stop at some value 
of 7]. Moreover, since a higher-order peak overtakes a lower-order one in height before the 
latter completely disappears, we indicate the highest peak in fig. |^ by the filled symbols. 

Referring to fig. ^ we can make the following remarks: according to our previous defini- 
tion, the fluid interacting by means of the potential given by eq. (|I|) is normal for packing 
fractions 7/ < 7/* = 0.50. Indeed, in this regime on the one hand the first maxima follow 
the scaling ki oc p^^^ and on the other hand the higher-order maxima kn are located in 
positions kn = nki, both features being a manifestation of the existence of a single length 
scale p~^^^ in the structure of the system. Above the overlap density p* the scaling breaks 
down. However, if we extrapolate the straight line with slope 1/3 which characterises the 
normal regime to higher densities, we find that it passes precisely through the main maxima 
at the magic packing fractions rim- Indeed, at those values of the density the length scales 
p~^^^ and a are commensurate and we have an accidental scaling of the main peak of S{k) 
with the 1/3-power of the density. However, the higher-order peaks of S{k) are still not 
located at integer multiples of the first one and we are dealing with an unusual fluid for all 
densities exceeding the overlap density p*. 




77i = 3, 4, 5, ■ ■ ■ 



(3) 



6 



IV. HARD SPHERE MAPPING AND THE FREEZING TRANSITION 



The structure factors obtained for the whole range of densities have the characteristic 
that the maximum height of the main peak in S{k) never exceeds the value 2.8. According 
to the empirical Hansen- Verlet criterion [jl2| , a liquid freezes when the first maximum of the 
structure factor reaches a value 2.85. It was subsequently shown that indeed freezing sets in 
if the order in the liquid phase, as measured by the first maximum of S{k) exceeds this quasi- 
universal value [^,0. Hence, we have a first indication that for / = 32 the system does not 
freeze. Naturally, the same must also hold for smaller /-values as the logarithmic core is then 
even softer. The absence of freezing was also confirmed in our numerical simulations, where 
for all values of r] which were simulated [r] < 12.00) the system remained in a liquid-like 
configuration, with any order parameter of a hypothetical crystalline solid vanishing. 

Here, we want to apply a simple Hard-Sphere mapping procedure in order to corroborate 
the above results for / = 32 on the one hand and establish approximately the 'critical arm 
number' fc on the other, which is defined as follows: for f < fc the system remains fiuid 
at all densities, but for f > fc there is at least one region in the density domain where the 
system is crystalline. 

Following an idea of Kang et al. [|l5l, we define an effective Hard-Sphere diameter ansi 



crudely representing the particle repulsion embodied in the pair potential of eq. (|I|), as 
follows. The pair potential V{r) is divided into a short-ranged reference potential Vo(r) and 
a longer-ranged perturbation potential W{r) at a suitably chosen break point A. Thus, we 
write 

V{r) = Vo{r) + W{r), (4) 

where Vo(r) is given by 

Vq^t) = e(A - r) [V{r) - F{r)] . (5) 

Here, 6(r) is the Heaviside step function, and F{r) will be defined soon. According to eqs. 
(I) and (H), the perturbation potential W{r) is given by: 

W{r) = e{X-r)F{r) + e{r-X)V{r). (6) 

Having splitted the potential in this way, it is now possible to calculate the free energy of the 
system by a simple Hard-Sphere perturbation theory [r^,|TB|. In this theory, the potential 



Vo(r) is used to calculate an effective Hard-Sphere diameter ans, while the longer-ranged 
potential W{r) is treated in a mean- field approach. In the following, we are not going to 
trace out calculations of free energies, but only to calculate ans, which is density-dependent 
when the specific choice for A and F{r), applied successfully in Refs. [|T3],|TB|, is used. This 
density-dependence arises from the identification of A with the nearest-neighbour distance 
a fee of the fee structure, i.e., 

A = «/cc-(^yj . (7) 
Furthermore, the function F(r) is chosen as 
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F(r) = V^(A) - 



dV{r) 
dr 



(A - r) , (8) 



guaranteeing that both Vo(r) and dVQ{r)/dr are vanishing at r = A. Having specified A 
and F{r), ans is then calculated from Vo(r) by the well-known Barker-Henderson (BH) 
approximation p]7| 

CTHS = dr 1- e-^^°(^) , (9) 



or, alternatively, by a scheme proposed by Weeks, Chandler and Anderson (WCA) ||T5|,|T6|,|T8 



We will only present results from the BH approximation, since the corresponding WCA 
results are practically identical to the latter, thus leading to the same conclusions. 

After obtaining cths from the BH or WCA scheme, we calculated the effective Hard- 
Sphere packing fraction rjns = ^P^'hs^ which is in general density-dependent. We can hence 
plot rjHS as a function of the 'true' packing fraction 77, depending furthermore on the value of 
/ as additional parameter. In fig. ^, we show results for rjns for five different functionalities 
/. Obviously, rjus is increasing linearly with rj for 77 ^0.1, and it reaches a maximum value 
depending on / at certain packing fractions 0.4 ^77 ^0.7, which are again depending on /. 
For rj > 0.74, rjHs remains constant for all densities. 

We will now briefly explain this behaviour ofT]Hs{v)y before switching to some conclusions 
that can be drawn from flg. ^ concerning the freezing transition of star polymers. For small 
rj, where the break point A is located at distances larger than the range of the pair potential, 
eqs. @, (IP and @ lead to a density-independent (Jhs, since the pair potential V{r) is 
density- independent. Hence, rjns scales linearly with 77 in this regime. However, when A 
reaches distances where the potential is remarkably different from zero, this linear scaling 
is no longer valid, since, according to eqs. and (|), cths is now a decreasing function 
with increasing density. This fact materialises in a decreasing slope of the function 7/^^5(77), 
leading even to a existence of a maximum in rjfjsiv)- Having these arguments in mind, the 
surprising fact that rjns attains a constant value for 7/ > rjc = 0.74 can be explained as 
follows. Exactly at the cross-over packing fraction 7/c, the break point A is located at the 
corona diameter a and thereafter, for r] > rjc "we have X < a. Therefore, for rj > t]c the 
Yukawa part of the interaction potential V{r) is irrelevant for the calculation of cths, as the 
reference potential Vo(r) is purely logarithmic. Using eqs. (|l|) and (P), it follows that in this 
regime ans scales linearly with A, i.e., 

(THS = A{f) A, (10) 

with the constant A{f) depending only on the arm number /, but neither on a nor on A. 
Since A scales with p~^^^, eq. (0) leads directly to the observed density-independence of 

VHsiv) for V > Vc [HI- 

The logarithmic pair interaction at hand is the only one showing this feature. Indeed, 
in order to have ans oc A, the integrand in eq. (^ must be a function of r/A only, i.e., 
the length scale a must drop out of the expression for Vo(r). Let us assume that the pair 
potential V{r) is given by some function R{r/cr] /) for r < A. If the reference potential Vo(r) 
has to depend on r/A only, the function R{x; /) must satisfy the following relations, as is 
clear from eqs. (1^) and (W) above: 
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R{r/a;f)~R{X/a-f) = R{r/X;f) and dMr/a^^l^ 

ar r 

where R{x; /) is some other function. The above conditions (which are, in fact, not 
independent but equivalent to each other) are fulfilled only by the family of functions 
R{x; /) = C{f) \n{x) + D{f), with C{f) and D{f) arbitrary. 

We now turn to the main reason for our interest in the values of the effective Hard-Sphere 
packing fraction rjns- Since the mapping onto an effective Hard-Sphere system gives quite 



reliable results in predicting the freezing transition |15,16|, we use the value of rjHs as an 
indication for the existence of a freezing transition in the original system. Hard Spheres 
freeze in a fee structure at rj^s'^ = 0.55 [^], and we therefore take this specific value to 



explore a possible freezing transition of star polymers from fig. ^. Obviously, for / < 32, 
rjHS never exceeds 0.55 for all rj, leading to the conclusion that the system remains fluid at 
all densities. For / > 32 on the other hand, rins attains values larger than 0.55 at least 
in a limited window of packing fractions rj. Consequently, there is a critical arm number 
fc ~ 32, meaning that the system never freezes for / < fc, but freezes presumably at least in 
a limited range of densities for / > fc- Our previously described finding that star polymers 
with / = 32 did not freeze in all our computer simulations is therefore consistent with this 
crude Hard-Sphere mapping procedure. 

As can be further seen from fig. |^, there is a range of arm numbers fc < /^64, where 
rjHs exceeds 0.55 only for 0.2^77^0.7, which implies a 'reentrant-melting'-phenomenon at 
T] ~ 0.7, i.e., a transition form a solid phase to a liquid phase if rj is increased above 
0.7. Here, it is worth mentioning that reentrant melting for star polymer solutions at high 
concentrations was predicted by Witten et al. on the basis of arguments arising from scaling 
theory [^. This prediction is independently verified here. Moreover, within our crude 
model, star polymers with /^64 show a liquid phase for 7^^0.2, followed by a solid phase 
for all higher densities. 



V. CONNECTION WITH EXPERIMENTS 

From the experimental point of view, the extreme values of the packing fraction that we 
have considered in section 3 are unattainable. At most, one can expect to observe a change 
in the behaviour of the first peak of the structure factor as the overlap polymer concentration 
(which corresponds to our overlap density) is crossed. Our prediction is that the height of 
the first peak of the star-star structure factor of a star polymer solution in a good solvent 
is not monotonically increasing with polymer concentration but rather it saturates at the 
overlap density and decreases thereafter. Based on general scaling properties of polymers in 



good solvents, Witten et al. predicted, more than ten years ago ||2T| , precisely that the peak 
of S{k) is largest when the separation between stars is of the order of the star radius. Here, 
we have confirmed this prediction quantitatively by employing a colloidal description of star 
polymer solutions. We expect that such effects should be visible in SANS experiments of 
dense star polymer solutions. To this end, stars with a labeled core should be used and 
the extent of the labeled part should be made as small as possible. This way, the effects of 
increasing concentration on the structure factor will not be masked by the form factor of 
the single star. 
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Additionally, it is natural to ask whether this nonmonotonic behaviour in the peak of 
S{k) is peculiar to the logarithmic form of the interaction inside the core or it can be seen 
for other functional forms of the pair potential as well. To test this, we have taken a toy 
potential which has a soft core for separations smaller than some length a and a crossover 
to a different functional form for larger separations. This toy potential reads as: 

0(r) _ [ A(r - a)V(ra) (r < a); ^^2) 



ksT [ (r > a), 

where A is some numerical constant which we can vary and controls the steepness of the 
interaction inside the core. We take now A = 10. Defining the packing fraction r] as in eq. (0) 
above, we have solved the HNC closure for a few different 77-values. The results are shown in 
fig. ^, where it can be seen that the same nonmonotonic behaviour of the main peak of S{k) 
is observed. Thus, the phenomenon is rather general and it relies on the existence of a soft 
enough core in the interaction. Such interactions are not uncommon in soft-matter physics. 
Hence, it would be of great interest if such anomalies could be experimentally observed. 

We finally comment on similarities of the logarithmic potential investigated in this work 
to a model introduced by Uhlenbeck and Ford in the early days of liquid-state theory ||22|| . 
This special model is defined by a Gaussian Mayer-function /(r) = 1 — exp{—f3V{r)) = 
exp(— ar^), corresonding to the pair potential l3V{r) = — ln(l — exp(— ar^)). As shown in 
Ref. all interparticle correlations can be calculated analytically within the framework 



of graph theory. For small distances r, the above potential potential reduces to (3V{r) ~ 
— 21n(Y^r), corresponding to the potential of eq. (|l]) with an arm number / = (36/5)^/^ ~ 
3.7. Since this arm number is significantly smaller than /c, we expect the model system 
of Uhlenbeck and Ford to be liquid for all densities, corresponding to an equation of state 
analytical in the density. 



VI. SUMMARY AND CONCLUDING REMARKS 

Employing a pair potential which has been shown to describe correctly the SANS data 
of star polymer solutions in good solvents [Q, we made quantitative predictions regarding 
the behaviour of the structure factor as a function of increasing polymer concentration. In 
particular, we found that above the overlap concentration, star polymer solutions display 
features which are unknown for normal liquids, namely a breakdown of the p^/^-scaling of 
the positions of the peaks of S{k) as well as an anomalous evolution of the heights of these 
peaks: the lower-order peaks diminish in height and the higher-order ones grow. 

Furthermore, we applied a Hard-Sphere mapping and used it as a crude diagnostic tool 
in order to make preliminary investigations on the topology of the phase diagram of star 
polymer solutions. The most striking feature of the Hard-Sphere mapping procedure is, 
in our view, the fact that the effective HS packing fraction rjus remains constant at high 
densities and depends only on the functionality /. We have shown that this characteristic 
is particular to the logarithmic interaction. We speculate that interactions which are softer 
than the logarithm will lead to an rjHs which will be a decreasing function of the density. 
This has relevance, in particular, to 'bounded' interactions, i.e., interactions which do not 



diverge at the origin, such as the Gaussian potential of Stillinger ||2^ or a simple model of 
penetrable spheres introduced recently by us p4| . 
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Another question of great interest is the phase diagram of star polymers. Detailed 
calculations, based on extensive computer simulations as well as perturbation theory are 
currently under way. These calculations allow for the evaluation of the Helmholtz free 
energies of various candidate crystalline phases and, subsequently, the comparison of the 
latter with that of the fluid phase and the construction of the phase diagram. Preliminary 
results are in full agreement with those presented here as far as the critical arm number fc 
is concerned; at the same time, they reveal a rich topology of the phase diagram as well 
as a variety of unusual crystal structures. The presentation of the phase diagram of star 
polymers will be the subject of a future publication. 
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FIGURES 

FIG. 1. The pair potential given by Eq. (|l|) for / = 18, 32, 64, 128, and 256 function of the 
centre-to-centre separation r. 



FIG. 2. The radial distribution function g{r) as obtained from simulations and from the RY 
closure for different r/-values. (a) ry = 0.10; (b) r] = 0.60; (c) r] = 1.50. 



FIG. 3. (a) The radial distribution function g{r) and (b) the anomalous structure factor S{k) 
above the overlap density. 



FIG. 4. (a) The radial distribution function g{r) and (b) the structure factor S{k) as obtained 
from simulations and from the RY closure at r] = r]2 = 3.40. For comparison, also the structure 
factor at r] = 2.00 is shown. 



FIG. 5. (a) The structure factor ^(A;) and (b) the radial distribution function g{r) for values 
of the packing fraction in the interval f?2 < ^ < %• 



FIG. 6. The structure factor S{k) at the magic values of the packing fraction r]m, for m = 2, 
3, 4 and 5. 



FIG. 7. Log-log plot of the positions of the various peaks of S{k) against the packing fraction 
ij. The filled symbols indicate the highest peak (see the text). The arrows indicate the locations 
of the magic packing fractions ijm and the dashed straight line has slope 1/3. 



FIG. 8. The effective Hard-Sphere packing fraction ijhs obtained by the procedure described 
in section 4 as a function of the true packing fraction rj for five different values of the functionality 



FIG. 9. The structure factor of a system interacting by means of the potential given by eq. (12) 
with A = 10 at various packing fractions r], obtained in the HNC closure. Notice the nonmonotonic 
behaviour of the maximum of S(k) with increasing density. 
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Watzlawek, Loewen and Likos, fig. 1 
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Watzlawek, Loewen and Likos, fig. 2(b) 




Watzlawek, Loewen and Likos, fig. 2(c) 




Watzlawek, Loewen and Likos, fig. 3(a) 




Watzlawek, Loewen and Likos, fig. 3(b) 




Watzlawek, Loewen and Likos, fig. 4(a) 




Watzlawek, Loewen and Likos, fig. 4(b) 
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